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^ We determine some basic properties of stars that produce spectacular gamma- 

CN ray bursts at the end of their Ufe. We assume that accretion of the outer portion 

^ of the stellar core by a central black hole fuels the prompt emission, and that 
Oh 

Q fall-back and accretion of the stellar envelope later produces the plateau in the 

c/3 X-ray light curve seen in some bursts. Using x-ray data for three bursts we 

a 

estimate the radius of the stellar core to be ~ (1 — 3) x 10^° cm, and that of the 

J>" stellar envelope tobe~(l — 2)xl0^^ cm. The density profile in the envelope is 

^ fairly shallow, with p ~ r^. The rotation speeds of the core and envelope are 

O 

1^ ~ 0.05 and ~ 0.2 of the local Keplerian speed, respectively. 

O 
00 

O Observations of gamma-ray bursts (GRBs) suggest that the activity at the center of these 



explosions lasts for several hours (1,2). The most compelling evidence is provided by three 



^ bursts (3) - GRBs 060413, 060607 A and 070110 - which show a sudden decline in their x- 

ray light curves (LCs) a few hours after the prompt burst (Fig. 1). The flux decline is by a 
factor of 10 or more and is much too sharp for the radiation to originate in an external forward 
shock (4); the most likely explanation is continued activity at the center of the explosion until 
at least the time of the decline. Additional evidence for continued activity of the central engine 
is provided by the x-ray flares seen in many GRBs (5-7), and also by those bursts whose x-ray 
& optical afterglow lightcurves (LCs) are mutually incompatible with a common origin (8, 9). 
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Figure 1: The top left panel (from Vaughan et al. 2006) shows the 0.3-10 keV band x-ray lightcurve 
(LC) of a typical long duration burst, GRB 050315. Four distinct phases are seen in the LC: the prompt 
GRB lasting for about a minute, an early steep decline lasting for ~ 10 minutes, a plateau lasting until 
about 10^ s, and a "normal" post-plateau decay. The top right panel (from Liang et al. 2007) shows 
x-ray and optical LCs for GRB 070110. The x-ray LC falls very sharply at the end of the plateau 
(t ~ 3 X 10'^ s), before switching to a normal decline, whereas the optical LC is a single power-law 
for the entire duration. The two lower panels (also from Liang et al. 2007) show the x-ray LCs of 
GRB 060413 and GRB 060607A. Note the complex x-ray LCs of these bursts, with flares, breaks, and 
plateaus. In contrast, the optical LCs of GRBs are typically smooth and simple, such as those of GRB 
0701 10 and GRB 060607A, shown in the right panels in blue. 



In fact, central engine activity is implicated whenever the observed flux variability time scale, 
6t, is much smaller than the time elapsed since the onset of explosion. The reason is that for a 
relativistic external shock causality dictates that 5t>R/2cT'^, and the time it takes for photons 
to arrive at the observer from a shock front at radius R that is moving with Lorentz factor T is 
also ~ R/2cT'^, i.e. 6t/t ~ 1 for external shocks. 

In this paper, we adopt the coUapsar model of GRBs (10,1 1), in which the inner part of the 
progenitor star collapses to a rapidly- spinning black hole (BH) and the remaining gas from the 
star accretes on to the BH and produces an ultra-relativistic jet. With a few plausible assump- 
tions, we show that x-ray observations of the three bursts mentioned above (GRBs 060413, 
060607 A & 0701 10) may be inverted to infer the structure and rotation rate of their progenitor 
stars. These particular bursts were selected for this work because they have prominent x-ray 
plateaus which are almost certainly the result of central engine activity. Thus for these bursts, 
we have information on the power generation at the center covering an extended period of time, 
which enables us to determine the core and envelope structure of the progenitor star; only three 
bursts met this strict requirement of central engine activity dominated x-ray plateau. Some of 
the results we find regarding progenitor star properties - especially the core structure - apply to 
a much larger sample of bursts, as discussed at the end of the paper. We also note that optical 
emission was detected for two of the three bursts (Fig. 1). The optical LC is consistent with 
origin in the forward shock (FS). Moreover, for a reasonable set of parameters the amount of 
x-ray flux produced in the FS is found to be ~ 10^ times smaller than the flux observed during 
the plateau. Thus, there is consistency between a FS origin for the optical LC and central engine 
activity for the x-ray emission. 

The panel on the left in Fig. 2 shows a schematic ORB x-ray lightcurve, which is based on 
the data shown in Fig. 1. The panel on the right outlines the basic features of our model and 
shows how the four phases of the LC are connected to corresponding zones in the progenitor 



fall to the center, and it is approximately equal to the free fall time (12), tfb ~ 2{r'^ jGMy, , 



star. 

We assume that the star has a core-envelope structure, as is common in stellar models. The 
bulk of the mass is in the stellar core. 

Some of the mass in the core collapses directly to form a spinning BH and the rest accretes 
on the BH to produce the prompt GRB emission. Surrounding the inner zone is a second zone, 
which represents the transition region between the core and the envelope. This zone has a 
steeply falling density profile, and there is correspondingly a rapid decline in the x-ray flux. 
The outermost zone is a relatively low-density stellar envelope. Accretion of this gas produces 
the plateau in the LC. The drop in the x-ray flux at the end of the plateau, which can be sudden, 
corresponds to the time when the outermost layers of the envelope are accreted. 

The fall-back time is the time it takes a parcel of gas in the progenitor star at radius r to 

where M^ is the mass interior to r, and all times are in the frame of the GRB host galaxy. Let us 
assume that the accretion time tacc of the gas in the disk is smaller than tfb (see below). Then, 
at a given time t = 10^t2 s, the central engine will accrete gas that has fallen back from a radius 
r = 10^°rio cm, where 

no ~ l.5tf M^gi , (1) 

and Mbh,i is the BH mass in units of IOMq. We are assuming here that the interior mass Mr 
is dominated by the BH. Let us suppose that the gas at radius r in the progenitor star has an 
angular velocity equal to a fraction /^ of the local Keplerian velocity: 

nir) = Ur)n,ir) = Ur) {GMr/r^f/\ (2) 

Then, the specific angular momentum j* = 10^®j*,i8 cm^ s^^ of this gas is given by 

j;,i8(r) ^ 3.8 M^ginf ^(r) ^ 4.74/'M^gi^(r). (3) 
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Figure 2: The panel on the left shows a schematic x-ray lightcurve with the following four segments: a 
prompt emission phase, a steep decline phase, a plateau phase, and a post-plateau phase. For the three 
GRBs considered in this work, the last phase has a steep and sudden decline. The panel on the right out- 
lines our proposal for how the different segments in the LC are related to the accretion of corresponding 
zones in the progenitor star The radii (r) and spin parameters (Jq = rj/ilfc) of the various zones are 
estimated from the x-ray data. 



The gas will fall back toward the black hole until its angular momentum allows it to become 
centrifugally supported, at which point it becomes part of a thick accretion disk about the black 
hole. Its further in fall will then occur on the time scale at which viscous forces in the disk allow 
the orbit to decay. 

The viscous accretion time of the fall-back gas after it has circularized is approximately 
equal to tacc ~ 2/afifc(^d), where a = 0.1q;_i is the viscosity parameter, and r^ ~ J'^/GMbh is 
the radius at which gas with specific angular momentum j,,, circularizes. Here we have assumed 
that the accretion disk is geometrically thick, which is generally true at the high mass accretion 



rates expected in a coUapsar (13). We thus have 

4cc ~ 10 a:l/B(r)tft ~ 10 aZljl^sMsli s ■ (4) 

The first relation shows that the accretion time will be shorter than the fall-back time so long as 
/f7< 0.4 (for a -0.1). 

The jet power from the central engine is determined by the mass accretion rate Mrh on to 
the BH, which is usually a small fraction of the mass fall-back rate Mfb on to the accretion disk 
(the majority of the fall-back mass is expelled in a disk wind as it spirals in toward the BH [14]). 

If the mass fall-back rate decreases suddenly at some time to, e.g., because of the transition 
from the stellar core to the envelope or when gas at the outer edge of the envelope has fallen 
back, then it can be shown by the conservation of mass and angular momentum that the accretion 
rate on the BH will decline with time as (15) 

MBH(t)~MBH(to)[l + 1.5(t-to)/4cc]"', t>to. (5) 

The exponent outside the square brackets is model-dependent and varies between —4/3 and 
—8/3, depending on details; for definiteness, we choose —2. If tacc ^ ^o^ the effect of the 
above time dependence is that, in a log-log plot, the jet power will initially drop by a large 
factor ~ (to/^acc)^ within a time ~ to^ and the power will then transition to an asymptotic 
decline oc t"^. This feature in the predicted LC plays an important role in our model. It is how 
we explain the sudden decline in the x-ray LC at the end of the prompt GRB and at the end of 
the plateau. 

We apply these scaling relations to derive constraints on the stellar structure and rotation rate 
of GRB progenitors. A number of the features seen in GRBs can be explained naturally within 
the coUapsar model as corresponding to the accretion of different portions of the progenitor star 
onto the central BH. Here we quantify the size, density, and rotation rate of the stellar core and 
envelope, as inferred from the observed timescales discussed above. 
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Observations with the x-ray telescope aboard the Swift satellite show that the flux in the 
0.3-10 keV band of a typical long duration GRB is roughly constant for about 20 s (Fig. 1, top 
left panel), after which the x-ray LC undergoes a rapid decline as ~ t^^ or faster for about 5 
minutes. Assuming that it takes ~ 30 s for the BH to form, and adding to this the burst duration 
of 20 s, we estimate that the end of the main burst is approximately 50 s after the initiation 
of core collapse. By equation ([T|), the gas that falls back at this time comes from a radius 
~ 9 X lO^Mgjj ]^ cm. This radius must correspond to the outer edge of the innermost zone in the 
schematic stellar model shown in Fig. 2. 

The subsequent steeply declining phase after the initial burst goes from t ~ 50 — 300 s. 
In our model, this portion of the LC is associated with accretion of gas from the second zone 
in the star, the transition region between the core and the envelope where the density has a 
steep decline. Using the same argument as before, the outer radius of the second zone must be 
~ 3 X lO^^MBgi cm. 

A rapid decline in the x-ray LC requires a rapid decline in Mbh (faster than ~ t~^), which 
is possible only when the disk is able to adjust quickly to a decrease in the mass-fall-back rate 
Mfb. For an accretion flow to respond quickly to rapid changes in Mfb at time t we must have 
^acc -C t (see Eq. |5]). This condition provides an upper limit on the rotation speed in the core 
region of the progenitor star. The x-ray flux typically drops by a factor of ~ 10^ during the 
rapid decline phase at the end of the prompt burst (Fig. 1). Therefore, by Eq. |5} we require 
(250/tacc)^ ~ 10^ or tacc ^ 8 s at t ~ 50 s. From Eq. 4, this means we must have fn < 0.2a_i at 
r ^ 10^° cm. 

We can obtain a lower limit on fa in the stellar core by requiring that the angular momentum 
should be sufficiently large that the fall-back gas is able to form an accretion disk, i.e., rd > 3Rg, 
where Rg = GM^u/c^ (note that the radius of the innermost stable circular orbit is 2.3Rg for 
a BH with a spin parameter of 0.9). Since during core collapse a particle initially at r ends up 



at a radius rd ~ r[fQ(r)]'^, where it is centrifugally supported, the condition r^ > 3Rg yields 
fn{r) > 0.02. An additional constraint on the core rotation rate is provided by the total energy 
of ~ 10^^ erg observed in a typical GRB. Numerical simulations show that the efficiency for jet 
production from a rapidly rotating BH is about 1% (16), and therefore the total mass accreted 
by the BH should be about 0.5Mq. This suggests that much of the mass within r ^ 10^° cm 
collapses directly to the BH, and therefore fn{r) is not much larger than 0.02; a reasonable 
estimate is /n ~ 0.05 for r < 10^° cm. 

A plateau in the x-ray LC is seen typically from ~3xl0^ — 3xl0^s (time in host galaxy 
rest frame). Gas falling on the central accretion disk during this time interval arrives from 
3 X lO^^Mg^f^ cm < r < 1.5 X lO^Ufg^f^ cm (Eq. fl|). We suggest that the gas accreted during 
this phase comes from the envelope of the progenitor star, which means that the outer radius of 
the envelope must be ~ 1.5 x 10^^ cm. In principle, the plateau could be produced by supernova 
ejecta that failed to escape. However, such a model is unlikely to give the sharp cutoff which is 
seen in the x-ray LC at the end of the plateau in the three GRBs considered in this paper. For 
these GRBs at least, the plateau must arise from fall-back of the stellar envelope, which has a 
well-defined outer edge. 

The density profile in the stellar envelope can be determined from the rate of decline of the 
x-ray flux during the plateau: f^ oc t^^, (5 ~ 0.5. If the specific angular momentum of the gas 
in the envelope is constant, as might be the case if the envelope is convective, then the required 
density profile is p(r) oc r^^^^+'^^f'^ ~ ^-2.2 ^^^ ^ ^ g g (15)_ For a convective stellar envelope 
that is not very massive, we expect p oc r^'^/^ when gas pressure dominates and p oc r~^ when 
radiation pressure dominates. The density profile we infer, assuming a constant specific angular 
momentum envelope, lies between these two limits. If the specific angular momentum increases 
with increasing radius, then a shallower decline of density is indicated; for instance, p oc r^^^ if 
the specific angular momentum increases with distance as ~ r^^^. A solid-body rotation profile 
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for the envelope is ruled out because it requires an unphysical density structure in which p is 
almost independent of r. 

The x-ray flux at the end of the plateau in GRBs 060413, 060607 A & 070110 faUs off 
sharply by a factor ~ 10^. The fall-off occurs at t ~ 3 x 10'^ s and extends over a time 5t ~ t. 
This steep fall-off requires the accretion time to be sufficiently short, tacc ~ V^O ~ 300 s, which 
implies that jd^ig < Sai/fMggi and ^(r) < Q.2al^ at r ~ 1.5 x 10" (17, 18). 

A lower limit for /^ in the stellar envelope is obtained by the requirement that the in-falling 
gas should be able to form an accretion disk at the center; this implies fu > 0.01. Another 
constraint comes from the fact that, apart from the flares discussed below, the x-ray LC during 
the plateau is usually quite smooth. This suggests that tacc is fairly large (which plays the role 
of a smoothing time scale), and so fn is probably not much less than the upper limit of ~ 0.2. 

After the sharp fall-off at the end of the plateau, the LC is expected to decline as ~ t^ if 
the x-ray emission is dominated by jet luminosity associated with the debris disk (see Eq. [5]). 
However, at these low flux levels the observed x-rays might be dominated by emission from 
shock-heated circumburst gas, i.e. external shock emission. In this case the flux will decline as 
~ t-\ e.g., GRB 0701 10 (Fig. 1). 

For those GRBs where the x-ray LC makes a smooth transition at the end of the plateau 
phase, fn should be such that tacc ~ t. For these systems, we require /^ ~ 0.4 near the outer 
edge of the progenitor star. 

Rapid flares with short rise times are often seen during the plateau phase of the x-ray LC, 
and these provide additional constraints. A flare was seen at the beginning of the x-ray plateau 
in GRBs 060413 and 060607A, and at the end of the plateau in GRB 070110 (Fig. 1), with 
a rise time of order O.lt. Assuming the flares are produced by a disk instability, the rise time 
should be a factor of a few larger than the instability time scale tinst- For a viscous instability 
^inst ~ iacc whereas for a dynamical (e.g., gravitational) instabihty ti^st ~ ^k^ ~ a^acc- Taking 



1/3 

the flare rise time to be ~ Stjnst ~ 0. It we find that /c ~ 0. la_i in the steUar envelope if flares 
arise as a resuh of viscous instability, and /q ~ 0.3 for a dynamical instability origin for flares. 
The amplitude of the flare fxfiarc/ fx, plateau can be at most ~ tacc/tmst- Thus, even in the limit of 
a dynamical instability, the amplitude is limited to ~ a^^ ~ 10. Much larger flare amplitudes 
(e.g., in GRB 060526 [5]) might suggest an unusually small value of the viscosity parameter a. 
Alternatively, these flares may be caused by a sudden increase in the mass fall-back rate, though 
such an event is not easy to visualize in our model. 

The decay time of flares is expected to be of order tacc, as this is the time scale on which a 
transient enhancement of the accretion rate, regardless of its origin, will subside. A noticeable 
difference between the x-ray plateau and the prompt emission is that LCs are typically more 
variable during the burst, and this raises a question as to why the central engine behaves differ- 
ently during the two phases. Part of the difference is that instability timescales are longer when 
the outer envelope of the star, with a larger specific angular momentum, is accreted. Another 
factor is that, during the plateau phase, the jet propagates through an already evacuated cavity 
and is less prone to fluctuating baryon loading. 

An upper limit on the stellar radius i?* is provided by the requirement that the energy ex- 
pended as the relativistic jet makes its way out of the star not exceed the energy initially injected 
into the jet. Using results in (19, 20) we find /?* < 5 x 10^^ cm. Our radius estimates are con- 
sistent with this limit. 

A possible way to get a handle on the mass of the GRB progenitor star is via the total energy 
produced in the explosion. The mass accreted during the prompt burst was estimated in §2.2.2 to 
be ~ 0.5Mq. The energy release during the plateau is typically about 10% of that in the prompt 
burst, i.e., ~ 10^^ erg (assuming the same beaming factor and efficiency as during the initial 
burst). Therefore, the mass accreted by the BH is only ~ O.O3M0. Even after allowing for the 
fact that only a small fraction ~ [r^/ RgY^"^ ~ 10^^ of the total fall-back mass actually reaches 
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the BH, the rest being carried away in a disk wind, we estimate that the mass of fall-back gas in 
the plateau phase is no more than a few Mq. The mass of the BH can be constrained because 
much of the mass within ~ 10^° cm collapses to the BH; for an evolved star this mass is about 
5Mq. Thus we estimate the mass of the GRB progenitor star to be M* ~ IOMq, plus whatever 
mass is ejected in the accompanying supernova explosion, which we are unable to constrain. 

Figure 2 summarizes our primary results on the properties of GRB progenitor stars. The 
three bursts we considered in this paper are notable in that they have a steep fall-off of x-ray 
flux at the end of the plateau phase. These bursts provide the most detailed information on 
the properties of their progenitor stars. Many of our arguments and results, especially those in 
which we use the prompt burst and subsequent steep fall-off to infer the properties of the stellar 
core, should apply to any long duration GRB. Similarly, our discussion of the x-ray plateau can 
be applied to other bursts that have plateaus in their x-ray LCs which are produced by continuous 
accretion activity i.e., those GRBs that show a simple power-law decline of the optical LC and 
a plateau in the x-ray; these comprise roughly one third of the long GRBs observed by Swift. 
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